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Abstract. We present an analysis of the properties of the cluster of galaxies RXJ 1821.6-1-6827, or NEP 5281, at a 
redshift z=0.816±0.001. RXJ 1821.6-1-6827 was discovered during the optical identification of the X-ray sources in 
^ ■ the North Ecliptic Pole (NEP) region of the 7? 05^4 T All-Sky Survey and it is the highest redshift cluster of galaxies 

i-^ ' of the NEP survey. We have measured spectroscopic redshifts for twenty cluster galaxies using the Keck-I and the 

, Canada-France-Hawai'i (CFH) telescopes. The value for the cluster velocity dispersion is av = 775lit^3 km s~^. 

' The cluster was also observed by XMM-Newton. Both the optical and X-ray data are presented in this paper. The 

cluster has an unabsorbed X-ray flux in the 2-10 keV energy band of F2-wkeV = 1-241q'23 ^ 10"^"^ erg cm~^ s~^ 
and a K-corrected luminosity in the same band of I/2-iofceV ~ 6.32^0 73 x 10** h^Q^ erg s~^ (90% confidence level). 
The cluster X-ray bolometric luminosity is Lbol,x = 1.35to 21 ^ 10*'^ Ko erg s"^ {Lbol,x = 1.17to \l ^ 10*'^ 
JH , h^Q^ erg s~^ in the concordance cosmology). The data do not allow fitting both metal abundance and temperature 

at the same time. The abundance is unconstrained and can vary in the range 0.28-1.42 Zq while the best fit 
X-ray temperature is T=4.71q'7 keV. This emission weighted X-ray temperature is a little lower, barely within 
the uncertainties, than the predicted temperature, T=6.34lg'35 keV, from the Lx — Tx relation of local clusters 
published in the literature. The optically measured velocity dispersion is consistent with the velocity dispersion 
expected from the ay — Tx relationship. We also examine the point X-ray source RXJ 1821.9-1-6818, or NEP 5330, 
located to the south east of the cluster which was identified as a QSO at z=1.692±0.008 in the ROSAT NEP 
survey. The X-ray source is well fitted by an absorbed power law model with Nh ~ 4.1I3 J x lO'^" atoms cm^^ 
and a photon index F = 1.67±0.12 typical of an active galactic nucleus. 

Key words, galaxies: clusters - general - individual: RXJ 1821.6-1-6827; X-rays: general - individual: 
RXJ 1821.9-f6818 



1. Introduction 

Clusters of galaxies are the most massive, collapsed structures in the Universe. They are the highest pea ks in a cosmic 
terrain driven by gravitational clustering and represent the manifestations of cosmic structure building llPeebleslll99i 
|Peacocklll999ll . The study of the properties of clusters at intermediate and high redshift, both individually and as 
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a population, is important since it allows different cosmological and structure formation models to be tested and 
constrained. The internal mix of components within clusters, as well as the space density of the most distant and 
massive clusters and the temperature distribution function of X-ray clusters, can be used to determine fundamental 
cosmological parameter s llHenrv fc Arnaudll99ll:IOukbir fc Blanchardll992lll997HEke et alJl99aHBorEfani et al.ll99fll 
l200lUHenrvir2 000. 2004: Vik hlinin et alJl2003(l . Through the detection of X -ray emission from the hot intracluster gas. 
X-ray cluster surveys are unbiased in the sense that they exclusively select gravitationally bound objects and are 
essentially unaffected by projection effects. There are only a few k nown clusters of galaxies beyond redshift unity, 
most of them were found in X-ray survey samples (see iRosati et al.l l2004) . Even if we are far from having very large 
complete samples of distant (z > 0.8) clusters, our capabilities to discover and observe distant objects are rapidly 
increasing as more accurate and sensitive observing techniques become available. 

In recent years we have been involved in the study of the deepest region of the ROSAT All-Sky Survey, at the 
North Ecliptic Pole (fHe nrv et ai]l200ltl\^es et alJl200l|) . to produce a complete and unbiased X-ray-selected sample 
of clusters of galax ies. The resulting NEP cluster sample has been used to investigate the nature of cluster eyolutiq n 
ijGioia et alJ l200lj) and to explore the potential implications for large-scale structure models i Mulhs et al.ll20oi|) . 
RXJ 1821.6-1-6827, or NEP 5281, was detectd in the ROSAT^YjP survey as a 4.5cr source with only 4q±9 ne t photons 
in a vignetting corrected exposure of 5519s. The identification program of the survey ijGioia et alJl2003j) revealed 
the cluster nature of the source. Optical spectroscopic observations with the CFH and Keck-I telescopes placed the 
cluster at a redshift of z=0.816±0.001. RXJ 1716.6-^6708 is the most distant cluster of galaxies in the ROSAT NEP 
survey. Follow-up l^JslUA- Newton observations were also performed, and are reported here together with the optical 
observations. Note that this is the second z=0.81 cluster discovered in the NEP survey. The first one, RXJ 1716.6-1-6708 
does not appear in the final complete NEP sample because it did not meet the selection criterium 
for source count rate of signal-to-noise ratio > 4ct. 

Section 2 describes the optical and X-ray data acquisition and analysis. Section 3 details the results obtained and 
discuss these results. A brief summary is given in Section 4. Throughout this paper, we assume an Einstein-de Sitter 
model Ho — 50 km s~^ Mpc^^, VIm — ^, and 0,^ — for direct comparison to previous work in this field, but we also 
repeat the calculations in the current cosmological concordance model (/i, JIa/, f^A)=(0.7, 0.3, 0.7). At the redshift of 
the cluster, the luminosity distance is 5.62 h^Q^ Gpc, the angular size distance is 1.70 h^Q Gpc, and the scale is 8.26 h^p^ 
kpc per arcsec. Throughout the paper quoted uncertainties are 90% confidence levels for one interesting parameter. 



2. Observations and Data Analysis 

In this section we present the optical spectroscopy for the cluster galaxies performed at the CFHT and Keck-I and 
the X-ray follow-up observations of RXJ 1821.6+6827 acquired with XMM-Newton. We first describe the optical 
spectroscopy, and then present the X-ray data which show that NEP 5281 has a low temperature, lower than expected 
from its bolometric luminosity, but commensurate with its velocity dispersion. We also present the optical and X-ray 
data on the point source to the south cast of the cluster, namely RXJ 1821.9+6818, or NEP 5330, which is identified 
with a QSO at z=1.692. 



2.1. CFHT and Keck-I Spectroscopy 

RXJ 1821.6+6827 was observed on July 3, 1997, at the CFHT with the MOS instrument using the STIS2 2048^ and 
the O300 grism in multislit mode. The wavelength coverage was from 4000 to 10000 A and the pixel size was 4.8 A 
pix~^. The slitlet wid th of 1.5" provided a s pectral resolution of 17 A. The data were reduced using the MULTIRED 
package developed bv lLe Fevre et al.l l)l995|) . Seven galaxies from the CFHT data were identified as cluster members 
from the Call break only (one has [Oil] in emission). NEP 5281 was also obs erved at the Kcck-I telescope on June 
22, 2001, with the Low Resolution and Imaging Spectrograph fLRIS. Ioke et al. 1995') in slit-mask mode. With both 
spectrographs in operation, LRIS_R and LRIS_B, it is possible to obtain multi-object spectra covering the entire optical 
window in one integration. However we used only data from the red arm given the redshift of the cluster. The 600 
1/mm grating blazed at 7850 A, to target [Oil] plus Ca II break and CN at the redshift of the cluster. Combined with 
our grating angle the wavelength coverage was approximately 6500 A to 9500 A, and the pixel size was 1.25 A. The 
GG495 filter was used to eliminate the overlapping second order spectrum to avoid any contamination blueward of 
9500 A. Given the fact that galaxies are extended objects, even at this redshift, we used a slit width of 1.4", which 
gives a spectral resolution of roughly 6.5 A. For the selection of the objects to spectroscopically observe, we used deep 
two color (B and I) images previously taken by us at the University of Hawai'i (UH) 2.2m telescope. The UCSCLRIS 
package designed by Drew Phillips and collaborators at Lick Observatory was used to prepare the slit-mask files. We 
designed two masks at different position angles in order to cover the central part of the cluster. Only data from the 
first slitmask were used since the second slitmask exposure was affected by heavy cirrus and bad weather. Thus we 
do not have a complete coverage for the galaxies in the central region of the cluster. The first mask was exposed for 
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Fig. 1. Distribution of velocities for the 20 galaxies showing a low velocity tail at A„ —6000 km s"^. The curve 
shows a Gaussian distribution centered a.t < v >— 244, 736 km s~^ with width given by the derived ay — 1407 km 
s~^ (not corrected to the first order for cosmological redshift) and normalized for the 18 galaxies considered as cluster 
members. 



10,800s. After the mask exposure, flat fields and arc calibrations exposures were also obtained. The data have been 
reduced using the standard IRAF package routines for 2-D spectra images. The two-dimensional spectra had to be 
straightned using the software package WMKOLRIS kindly provided by Greg Wirth at the W. M. Keck Observatory. 

Fifteen galaxies out of 26 observed objects were recognized as cluster members on the basis of Call H and K, CN 
and, in some cases, the G band absorption features. Only one of the 15 galaxies observed had [Oil] in emission. Two 
of the seven galaxies observed at the CFHT are common to the Keck-I observations. In the end, 15 unique objects 
observed by Keck plus 5 unique objects observed by CFHT turned out to be cluster member galaxies. These twenty 
galaxies are listed in Table 1. For each galaxy the (J2000) coordinates, measured velocity, la error, and redshift are 
given. In the last column the main spectral features are also noted. 

The velocity histogram for the 20 member galaxies is shown in Fig. ^ There is a low-veloc ity extension in the 
histogram at 238,500 km (2 galaxies) that our 3ct iterative clipping algorithm (following iDanese et aLlFlOSOl) 
excludes from the computation of the cluster velocity. From 18 accepted cluster members, and taking into account 
the errors on the redshift, we obtain a mean velocity < v >— 244, 736 ± 325 km s^^, an average redshift < z > = 
0.8163 ± 0.0011, and a dispersion along the line of sight cry = 775li[^3 km s^^ (corrected to first order for cosmological 
redshift by dividing by 1-l-z). Without the 3a clipping the mean velocity for the 20 galaxies is < >= 244, 249 ± 384 
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Fig. 2. Figure 2 can be viewed and downloaded from this http URL'http: / / www.ira.cnr.it'^gioia/PUB / publications. htmlj 
The figure shows a 1800s I-band image of RXJ 1821.6+6827 taken at the UH 2.2m telescope. The combined MOSl 
and M0S2 image was smoothed with a Gaussian to have 16 pixels per cell (corresponding to 8"). The energy range is 
0.3-8 keV band. Units for the smoothed overlayed X-ray contours are (6.5, 7.8, 9.1, 13, 18.2, 24.7, 32.5) 10~^ counts 
arcmin^^ over the background which correspond to (2.5, 3, 3.5, 5, 7, 9.5, 12.5)(T over the background. The circles 
indicate 14 out of 20 cluster members. The cD galaxy is marked as galaxy # 7. The other 6 galaxies at the redshift 
of the cluster fall outside the field of view (see Tabic 1). The image measures 1092 x 1092 pixels covering a field of 
4' X 4' {l.98h^Q X 1.98h^Q Mpc at z 0.8163). 



data and folded model 




channel energy (keV) 

Fig. 3. Binned X-ray spectrum and residuals for RXJ 1821.6+6827 . The spectrum was binned in such a way that each 
resulting energy channel has a signal-to-noise ratio of at least 4. The solid line represents the best fit MEKAL model. 
The top grey spectrum is the pn spectrum, while the lower and darker spectra are the MOSl and M0S2 spectra. 
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^, the average redshift is < 2 >= 0.8148 ± 0.0016, and the dispersion along the line of sight is ay = 11421^55 
1 



Using the ROSTAT software by Tim Beers and collaborators ijBeers et al 



obtain very similar results to those of the 3cr iterative clip ping algorithm by 



|l990l) Csee alsolBird fc Beer£lll993|) we 
Danese et alJ l|l980l ). Below the values 



obtained by the biweight estimator of scale, as suggested bv lBeers et al ] ili99nir for this number of galaxies, are given. 
The errors are the 90% confidence intervals. The average redshift is < z >= 0.8156 ± 0.0024 accounting for errors on z 
(< z >— 0.8147±0.0044 with no errors), and a (Ty = 577±231 km s~^, corrected to first order for cosmological redshift 
{(TV = 1048 ± 421 km s^^ corrected to first order cosmological redshift but without taking into account the errors 
on redshift), which are consistent with cither the 20 cluster member and the 18 cluster member velocity dispersions 
obtained with the previously adopted 3cr clipping algorithm. 



2.2. X/WM- Newton observations 

XMM-Newton l* Jansen et al.l200lh observed RXJ 1821.6+6827 as part of the GO program in three epochs, on November 
13 and 23, 2002, and on December 19, 2002, for a nominal exposure of 63 Ks with the European Photon Imaging 
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Camera (EPIC) pn lIStrMer et alJl200lh and for 73.5 ks with the EPIC MOS CCD arrays llXurner et al.ll2n0lh . The 
pn was operating in extended-full-frame mode with thin filter while the MOS was in the full-frame mode also with the 
thin filter applied. The pileup is not a problem given the low count rate of the X-ray source. 

Unfortunately the observations suffer from periods of very high background. Event files produced by the standard 
pipeline processing have been examined and filtered to remove the high background time intervals (using the version 
5.4.1 of the Science Analysis Software, SAS, and the latest calibration files released by the EPIC team). Only events 
corresponding to pattern 0-12 for MOS and pattern 0-4 for pn have been used^. Good time intervals were selected by 
applying thresholds of 0.35 counts s~^ in the MOS and 1 counts s~^ in the pn to the photons at energies greater than 
10 keV. At these higher energies counts come mostly from background. The net exposure times, after data cleaning, 
are 21.3 Ks and 22.0 Ks for the MOS 1 and MOS 2, respectively, and 10.6 Ks for the pn. For subsequent analysis 
background counts have been accumulated using nearby source-free circular regions. 

The three images, one for each instrument, taken on different epochs were summed using the SAS task merge. 
Given the position of the source at the center of the field, the different orientations of the instrument with respect to 
the sky coordinates should not affect the response of the instruments. Thus resulting response matrices and auxiliary 
files are the average of the three exposures, weighted by the different exposure times. 

Response matrices (that include the correction for the effective area) have been generated using the SAS tasks 
arfgen and rmfgen. The X-ray fluxes reported below are computed using the MOSl detector and calibration(s). From 
our data the normalization of the M0S2 is 5% lower than that of the MOSl whereas the normalization of the pn is 
24% lower. This apparently discrepant normalization value has to be attributed to the fact that the cluster falls on 
the CCD gap of the pn. 



2.3. Overall Cluster Properties 

Fig-Elshows the XMM-Newton smoothed contours of RXJ 1821.6+6827 overlayed onto an optical 1800s I band image 
taken at the UH 2.2m telescope. Both MOSl and M0S2 data were used to produce the X-ray contours in the 0.3-8 keV 
band. The data were smoothed with a Gaussian using cells of 16 pixels corresponding to a resolution of 8". The X-ray 
cluster is extended along the NE-SW direction. There are not enough X-ray counts to produce a more meaningful 
comparison between the optical and X-ray images. 

To measure the emission-weighted cluster temperature, net counts were extracted from a circle of radius 100" cen- 
tered at 02000 = 18'' 21™ 32.5'', 62000 = +68° 27' 45.0" which corresponds to the barycenter of the X-ray emission 
and is very close to the position of galaxy 7 (see Table 1). At the redshift of the cluster such a radius corresponds to a 
linear size of 826 h^Q kpc. A larger radius including more extended emission counts could not be chosen given the high 
increment in background counts. There are 708±35 net counts in the (0.4-4.0 keV) energy range in the MOSl detector, 
671±36 in (0.5-4.3 keV) in the M0S2, and 751±40 net counts in (0.5-3.4 keV) in the pn detector. X-ray properties 
are listed given in Table 2. The data were fitted in XSPEC (version 11.3.0) with a single temperature MEKAL model 
llMewe et al.lll98 5|) modif i ed by Galactic absorption, where the ratio between the elements are fixed to the solar values 
as in I Anders fc Greve'ss3 l)l989j) ^. 

The relative normalization of the three instruments was left free to account for possible intercalibration differences 
or mismatches. The absorption and redshift were kept constant while the tempera ture and abunda. nce were left free 
to vary. The absorption was fixed at the Galactic value of 5.24 x 10^" atoms cm~^ ((Elvis et al .11 19941) and the redshift 
was fixed at z=0.8163. The best fit with a thermal MEKAL model gives a temperature T=4.71q y keV and a best fit 
abundance Z=0.77 Zq, which is unconstrained by these data due to a non-detection of the iron K line. The mctallicity 
can vary in the range 0.28-1.42 Zq. The binned X-ray spectrum and best fit folded model arc shown in Fig. O where 
bins were chosen so as to have each resulting energy channel significant at the signal-to-noisc ratio of at least 4. The 
total number of photons in each bin allows the use of the statistics for the fit. The fit has a reduced of 0.90 for 
102 degrees of freedom. If one fixes the abundance to Z=0.5 Zq then the temperature increases a little (T=4.9lg g 
kev) but it is consistent with the previously determined temperature within the errors. The reduced becomes 0.90 
for 103 degrees of freedom. 

The unabsorbed flux in the 2-10 keV energy band is i^2-iofeeV = l-24to'23 ^ 10~^^ erg cm~^ s~\ (see Table 
3). While the unabsorbed X-ray flux in 0.5-2.0 keV band is i^o.5-2.ofeey = l-35lo:2i X 10"" erg cm'^ s-\ to be 
compared with the flux obtained with the ROSAT PSPC in 0.5-2.0 keV of Fx.^rosat = 1-02 ± 0.22 x 10"" erg 
cm~^ s^^. The K-corrected luminosity in the 2-10 keV energy band is L2-iokeV = (^■'^'^-o rt ^ l^^^ 1^50^ SI'S while 

^ see the XMM-Newton Users' Handbook 

^http:/ /xmm. vilspa. esa. es/external/xmm_user-Support/documentation/uhb/XM M_ UHB.htm\ 

These values for the solar metallicity have recently been superseded by the new values of iGrevesse fc Sauvall (Il998l) . who 
used a 0.676 times lower Fe solar abundance. However, we prefer to report metallicities in the units of Anders & Grevesse 
abundances since most of the literature still refers to these old values. 
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Fig. 4. Long-slit spectrum taken at the UH 2.2m of the QSO at z=1.692 associated to RXJ 1821.9+6818. The total 
integration time was of 2400s. The dashed Hncs indicate the positions of the CIV, Clll], Mg emission lines at the QSO 
redshift. Wavelengths of atmospheric absorption bands are also indicated. 



the bolometric luminosity is Lbol,x 
with XMM-Newton is L0.5-2.0keV = 4.04 

of LxMOSAT 



l-35to:2i X lO**^ Kq erg s^^. The luminosity in the 0.5-2.0 keV obtained 



+0.47 

0.58 



X 10'*'^ erg s"\ comparable with the ROSAT PSPC luminosity 

(2.91 ± 0.63) X 10*4 

^50^ erg s ^. In the concordance cosmology the K-corrected luminosity in the 



2-10 keV energy band becomes L2-iokeV — 5.481q;74 x 10"^^ h^g^ erg s ^ and the bolometric luminosity becomes 



Lbol.x - 1.17t^:i^ X 10^' hfo' erg s" 

As can be seen from Fig. |2 there are point-like sources both to the north-east and the south-west which are 
included in the 100" radius even if it is unclear if they belong to the cluster. The data were re-analyzed excluding the 
point sources from the measure of the cluster X-ray temperature. Freezing as before the Galactic absorption to the 
neutral hydrogen column density N/f — 5.24 x 10^° atoms cm~^, the redshift to z=0.8163 and leaving the abundance 
free to vary, a value of T=4.9j;Q g keV is determined for the cluster gas once the point sources are excluded. The fit 
has a reduced of 0.91 for 85 degrees of freedom. This value for the temperature is perfectly consistent with the 
temperature previously obtained (T=4.7^q'7 keV) within the uncertainties. 



2.4. RXJ 1821.9+ 6818, or NEP5330: the QSO at z=1.692 in the XMM-Newton field 

There are at least 20 point sources in the XMM-New ton image. We d iscuss here only one point source, 
RXJ 1821.9+6818, which appears in the NEP survey catalog l|Gioia et alJl2003^ as NEP5330, and which is not visible 
in Fig. 2 since it is outside the field of view. The source was identified with an AGNl (or type 1 AGN, a classifica- 
tion which includes either Seyfert 1 and QSO objects) at redshift z=l. 692+0. 008. The object was spectroscopically 
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Fig. 5. Binned X-ray spectrum and residuals for RXJ 1821.9+6818, the QSO in the XMM-Newton field. The spectrum 
was binned in such a way that each resulting energy channel has a signal-to-noise ratio of at least 5. The solid line 
represents the best fit power law model. The top grey spectrum is the pn spectrum, the MOSl and M0S2 spectra are 
the two lower and darker spectra. 

observed at the UH 2.2m telescope as part of the identification program of the NEP survey sources on June 18, 1998. 
The telescope was equipped with the Wide Field Grism Spectrograph and the Tek2048 CCD. We used the 420 1 mm^^ 
red grism and a 300/i slit (2.4") which provided at the f/10 focus a pixel scale in spectroscopic mode of 3.6 A pix^^, 
a spectral resolution of ~ 19 A FWHM, and a wavelength coverage of approximately 3900A-9500A. 

The classification with an AGNl has been done on the basis of the equivalent width of the emission lines (Wa ^30 
A) and on the FWHM (^4000 km s^^) of the broad permitted emission lines, indicating a QSO type object (see 
ICioia et alJ 2003 for more details on the optical classification). An optical spectrum is shown in Fig. 01 where the 
dashed lines indicate the positions of the emission lines at the QSO redshift. Wavelenghts of atmosphric absorption 
bands are also indicated. 

To examine the X-ray source associated to the QSO, a circle of radius 25" was centered at Q!2ooo = 
18'' 21"* 58.3^ (52000 = +68° 18' 41.2". There are a total of 609±25 net counts in the (0.4-3.8 keV) energy range in 
the MOSl detector, 632±26 in (0.5-4.6 keV) in the M0S2, and 874±30 net counts in (0.5-6.2 keV) in the pn detector 
(see Table 2). The data were fitted in XSPEC with an absorbed power law model given the nature of the source. 
The relative normalization of the three instruments was left free. The X-ray flux reported below is computed using 
the MOSl detector and calibration. The normalization of the M0S2 is 2% lower than that of the MOSl whereas the 
normalization of the pn is 3% lower. The three instrument normalizations are in better agreement than in the analysis 
of the cluster source implying that the high normalization discrepancy (24%) found previously for the pn with respect 
to the two MOS detectors is most probably due to the fact that the cluster falls in one of the gap of the pn. 

Both the absorption and the photon index (F) of the power law model were free to vary. The best fit F is 1.67±0.12 
and the best fit absorption is N/f ~ 4.1 Ij'q x 10^" atoms cm^^, consistent with the Galactic value of = 5.24 x 10^° 
atoms cm^^. The binned X-ray spectrum and best fit folded model are shown in Fig.[5lwhere bins were chosen so as to 
have each resulting energy channel significant at the signal-to-noise ratio of at least 5. The fit has a reduced of 0.92 for 
61 degrees of freedom. The unabsorbed flux of the QSO in the 2-10 keV energy band is F2-iokeV = (3.66±0.39) x 10^^"^ 
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erg cm ^ s ^ (see Table 3). The K-corrected luminosity in the 2-10 keV energy band is L2-iokeV — 5-0lo.7 x 10"*^ h^Q 
erg s^^. In the concordance cosmology the K-corrected luminosity in the 2-10 keV energy band becomes i2-iofeey = 
(5.3;o;5) X 10^5 hfg^ erg s'^. Note that the luminosity of RXJ 1821.9+6818 as m easured with XMM-Newton in 2002 
is approximately 2.5 times higher than that measured with ROSAT in 1990 fsee lMuUis et alJ 2004 for revised fluxes 
and luminosities for NEP AGN). Such a difference is not surprising given the ti me variability of AGN. 

RXJ 1821.9+6818 was also observed with the Very Large Array at 1.5 GHz bv lKollgaard et all l)l994|) who mapped 
a 29.3 deg^ region surrounding the north ecliptic pole to support the deepest portion of the ROSAT all-sky soft X-ray 
survey RXJ 1821.9+6818 was detected by the the VLA-NEP survey with a flux density of 182 mJy at 1.5 GHz. 



3. Discussion and Conclusions 

The main result of this paper is the presentation of the optical and X-ray data, namely the spectroscopic data 
acquired through Keck-I, CFH and UH 2.2m telescopes, and the X-ray temperature obtained with the detectors 
onboard the XMM-Newton satellite. We will focus the following discussion on the implications of measured quantities 
such as the emission weighted temperature and velocity dispersion for RXJ 1821.6+6827. As with other X-ray selected, 
distant clusters RXJ 1821.6+6827 does not have a spherically symmetric morphology in X -rays or in optical. Its X-ray 
morphology is elongated and very similar to other ~0.8 clusters like RXJ 1716.6+6708 (|GioiaeE^[l99i), another 
cluster foun d in the course o f the NEP identification program but which did not meet the Aa selection criterium 
for sources ijGioia et alJl2003l) . or to the better known Medium Survey cluster MS 1054-0321 (|Gioia et al .1 120041) . 
Differently from RXJ 1716.6+6708 and M S 1054-0321 that both show high velocity d ispersions (cry = 1522^111 km s'^ 
and av = 1150 ±97 km s-\ resoectivelv: iGioia et alJll999Hvan Dokkum et allbOOOl) . RXJ 1821.6+6827 has a velocity 
dispersion of ay = '775~!i\i^ km s""'^, rather typical of normal relaxed clusters. We are also dealing with a cool cluster 
with a temperature (T =4.7 keV) a little low er than the temperature of RXJ 1716.6+6708 (which Chandra measured 
to be T=6.6±0.8 keV: IVikhlinin et alllioO^ and even lower th an the X-ray temperature of MS 1054-0321 (which 
XMM-Newton measured to be T=7.2t°;^ keV: lGioia eraLll2004|) . 

A large number of authors (see Table 5 in Girardi et al. 1996, or Table 2 in Wu et al. 1998 for an exhaustive 
list of papers on the subject) have attempted to determi ne the ay—Tx using different cluster samples in order to 
test the dyn amical p roperties o f clusters (s ee among oth e rs iGirardi et alJ 1998t Donahue et al.lll999t iTran et al]ll999l: 



iGioia et alj ri999: Ho rner et al1ll99 9: Ros ati et al.ll2002l: lYee k Eningsonll2003t iLubin et alJl2004|r This relationship 



is physically meaningful since both the velocity dispersion of the galaxies and the temperature of the intracluster 
medium provide a measure of the overall mass of the system. If we characterize the ay—Tx relationship as 

kT 

"--^ gas 

(where fj, is the mean molecular weight of the gas, mp is the proton mass and k is the Boltzmann's constant) and if we 
assume energy equipartition between the galaxies and the gas in the cluster (/3=1), we obtain a velocity dispersion for 
RXJ 1821.6+6827 consistent with the measured one. We derive ay = 848!*;g7^ km s~^ from the above ay—Tx relation 
vs the measured value of ay = 775^113 km s~^. 

[Cirardi et al. (199^ have derived a best fit relation betwee n the velocitv dispersion and the X-rav temperature, with 
more than 30% reduced scatter with respect to other workers ((Edge fc Stewart! 1 99 lULubin & Bahcalil99.^lBird et alJ 
Il995t IWu et al.lll998l ). iGirardi et al.l l|l996l) have taken into account distortions in the velocity fields, asphericity of 
the cluster or presence of substructures to derive their best flt relation: \og{ay) — (2.53±0.04) + (0.61±0.05) log(T). 
If we insert the temperature of RXJ 1821.6+6827 in the above relation, the resulting velocity dispersion value, ay — 
874^83^ km s~^, is consistent, even if on the high side, with the measured velocity dispersion. 

The temperature of RXJ 1821.6+6827 is a little lower than predicted from its X-ray luminosity for local Lx—Tx 
relation. Understanding the evolution of the Lx — Tx relation is important not only to understand the physics 
behind the formation of galaxy clusters, but also because it provides a link between observations of clusters and 
derivation of cosmological parameters. There is an extensive literature on the correlation between these two basic 



among other "David et al."l99^ 'Fabian et al."l994^. 'Mushotzk> 


' & Schardll997*. 


'Markevitchlll998l 


Sadat et alJll998l: 


Arnaud & Evrard 1999: Fairlev et al. 2000: Borgani et al. 200lHNovicki et al.ll2002: Vikhlinin et al.ll2002t 




2004(). The Ly — Ty relation has been well studied at low redshift. Characterizing 


the Ljsf-Tx relation as 



Lx = LJ—^)" 

investigators have found for a values ranging from 2.33 to 2.88, all consistent with in the erro rs. There is no conse nsus 
yet on the evolution of the Lx — Tx relation with redshift (see among others: iBorgani et al l l200L:.Novicki et al...2o'oi 
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IVikhlinin et alJ l2002t *Lu mb et al .1120041) p robably due to the lack of large samples of galaxy clusters at cosmologically 
significant redshift. Vikhli nin et al.l l)2002j) have detect ed for the first time an evolutionary trend in the hx-Tx relation 
using Chandra data for distant clusters. More recentlv lLumb et alJ l)2004|) have also found an evolutionary trend of the 
relation using XM M-Newton d ata, albeit using a sample of only 8 clusters in a restricted redshift range (0.45< z <0.62). 
On the other hand Isttori et a l. (2004) detect only hints of negative evolution in the hx-Tx using a larger sample 
of 28 clusters observed by C handra and in a larger redshift range (0.4< z <1.3). We assume the parametrization 
by lArnaud &: Evrardl l)l999(l who analyzed a sample of 24 low-z clusters with accurate ASCA temperature mea- 
surements and absence o f strong cooling flows. Comparing the bolometric luminosity of RXJ 1821.6+6827 with the 
lArnaud & Evrardl l|l999l) best fit relation, log{hBOL,x) = (2.88±0.15) log(T/6keV) + (45.06±0.03), one would expect 
for RXJ 1821.6+6827 a temperature of 6Mt°oll keV that is a little higher, even within the uncertainties, than the 
best-fit XyiM- Newton temperature, implying that this cluster may fit better in an evolving hx—Tx relation. Given 
the XMM-Newton results, RXJ 1821.6+6827 does not show extreme X-ray luminosity or temperature values. 

We can also estimate the mass of RXJ 1821.6+6827 using the X-ray temperature value. With the assumptions 
that the mean density in the virialized region is ~200 times the cr itical density at the redshift of the cluster and that 
th e cluster is isotherm al (Evr ard et al.lll995 iDonahue et al.lll998|) . we can use the scaling law method as illustrated 
in lArnaud fc Evrardl ((1999j). From the simulations of Evrard, Metzler & Navarro (1996) for the mass-temperature 
relation one can estimate the virial mass within a radius r2oo(r) — 1-85(yq|j^)^/^(1 + z)^'^/^h^^ Mpc by using the 
equation 

« (1.45 X lOi5ft-iM0)(l + z)-3/2(^L^)i-5 

IQkeV 

From the IQAlsll-Newton temperature the virial mass is approximately Mvir = (l-S^o 4) ^ ^^^^ ^50^ ^'^Q within 
'"200 = 1-04 h'^Q Mpc, where the uncertainties on the mass reflect the uncertainties on the temperature. 



4. Summary 

We have presented observations for the z=0.8163 galaxy cluster RXJ 1821.6+6827 performed with the instruments 
on board the XMM-Newton satellite, and with the CFHT and Keck-I telescopes. The main result of the paper is to 
present new optical and X-ray data for the most distant cluster in the NEP survey. Both the temperature and the 
metal abundance have been left as free parameters in the fitting of the X-ray data with a model. Thus freezing the 
hydrogen column density to the Galactic value of 5.24 x 10^" atoms cm~^ and the redshift to the measured z=0.8163, we 
obtain a best fit temperature with a thermal MEKAL model of T=^.7t.\''\ keV, while the abundance is imconstrained 
and can vary in the range 0.28-1.42 Zq. This X-ray temperature is a little lower than predictions from the cluster 
X-ray bolometric luminosity Lbol,x = l-35lE5:2i ^ '^'^^^Kq erg 8"^ (Lbol,x = 1.17^^X1 ^ ^'^^^h^Q erg s^^ in the 
concordance cosmology) from the Lx — Tx relation of local clusters published in the literature. 

From 18 accepted cluster members we obtain an average redshift < z >— 0.8163 ± 0.0011 and a dispersion along 
the line of sight cry = 775~!i\i^ km s~^. This measured value is consistent with the velocity dispersion expected from 
the ay — Tx relationship obtained from different cluster samples. From the M - T relation we estimate a virial mass 
for RXJ 1821.6+6827 within the r2oo — 1-04 h'^Q Mpc on the assumptions that the mean density in the virialized 
region is ^200 times the critical density at the redshift of the cluster and that the cluster is isothermal. The estimated 
virial mass has a value of M„r = l-Stoi x 10" ^50^ '^'hich shows that RXJ 1821.6+6827 is a massive cluster at 
high redshift. 

The point source to the south, RXJ 1821.9+6818, identified in the NEP survey as a QSO at z=1.692±0.008, is well 
fitted by an absorbed power law with N/f = 4.1^3 ^ x 10^" atoms cm^^ and a photon index F = 1.67±0.12 which is 
typical for this class of objects. 
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Table 1. Spectroscopic data for the galaxies in the NEP cluster RXJ 1821.6+6827 



ID 


RA (J2000) 


DEC (J2000) 


cz 


Acz 




z 


Spectral features 


# 


h m s 


o 


/ // 


km 


km 












1 


io Zi 


lo.D 


-rOo 




244367 


120 


0, 


.8151 


n-f-iv, 


Dana, 


TJX pxr 


2 


1 ft 91 


01 ^ 




9Q ^7 


246585 


240 


0, 


.8225 


xi-f-iv, 






3 


1 Q 01 
io Zi 






9Q "^1 


246016 


170 


0, 


.8206 


xi-f-iv, 






4 


1ft 01 
lo zl 


OU.O 




9Q 9Q 

Atj Zy 


245856 


180 


0. 


.8204 


n-f-iv, 


a VinnH 
\j uaiiLi, 




5 


1 ft 91 
lo zl 


97 Q 






246705 


480 


0, 


.8229 


n-f-iv, 


Vjr uaiiCi, 


PN 


6 


1 ft 01 
io zi 


oo.o 




9R AO 

ZO '4:U 


246885 


210 


0, 


.8235 








7 


1 ft 01 

io zi 


oZ.O 




97 


244847 


180 


0, 


.8167 


xi-f-iv, 


C hnnH 
VI uaiiLi, 




8 


1ft 01 

lo zl 


Q/1 O 




97 ^'^ 


242808 


90 


0, 


.8099 


H+K, 


CN 




9 


1 ft 01 

lo zl 


Of .Z 


-rOo 


97 fil 


244428 


389 


0. 


.8153 


H+K, 


G band, 


pi\r 


10 


1 Q 01 

lo zl 


OO.O 




97 ^0 


243317 


389 


0, 


.8161 


H+K, 


G band 




11 


18 21 


35.2 


-hDo 


97 Q/1 


243288 


60 


0, 


.8115 


H+K, 


CN 




12 


18 21 


36.6 


+68 


27 21 


244996 


120 


0, 


.8172 


H+K 






13 


18 21 


37.5 


+68 


27 07 


246765 


150 


0. 


.8231 


H+K 






14 


18 21 


37.9 


+68 


27 00 


244607 


60 


0. 


.8159 


H+K 






15 


18 21 


18.2 


+68 


27 00 


238102 


840 


0, 


.7942 


[OH] 






16 


18 21 


47.7 


+68 


27 06 


242508 


30 


0, 


.8089 


H+K 






17 


18 21 


51.1 


+68 


27 17 


241639 


540 


0, 


.8060 


Call Break 




18 


18 21 


57.7 


+68 


27 25 


243378 


540 


0. 


.8118 


Call Break 




19 


18 21 


50.8 


+68 


26 09 


244397 


840 


0, 


.8152 


Call Break, [Oil] 


20 


18 21 


47.8 


+68 


25 18 


243498 


90 


0, 


.8122 


H+K 







Table 2. X-ray properties of RXJ 1821.6+6827 (NEP 5281) and RXJ 1821.9+6818 (NEP 5330) 



RA (J2000) DEC (J2000) Instr R Net Cts T Z (Z©) 

h m s ° ' " " keV r 



N/f Region 
xlO^° at/cm^ 



18 21 32.5 +68 27 45.0 



18 21 58.3 



+68 18 41.2 



MOSl 


100 


708+35 


M0S2 


100 


671±36 


pn 


100 


751+40 


MOSl 


25 


609±25 


MOS2 


25 


632±26 


pn 


25 


874±30 



4.7^ 



+0.65 



0.77 



1.67+0.12 



frozen'' 



4.1 



+3.1 



cluster 



QSO 



Nif fixed at the Galactic value of 5.24 x 10 atoms cm 
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Table 3. Fluxes and luminosities for the galaxy cluster and the QSO 



RA (J2000) 


DEC (J2000) 


F2-10fceV 


z 


La 


-lOkeV 


Lbol 


Region 


h m s 


o / // 


10"" cgs 




10** 


^50 Cgs" 


10*^ h^t 


cgs 










10** 


iiro cgs'' 


10*« 


cgs 



18 21 32.5 


+68 27 45.0 


^•^^-0.23 


0.816+0.001 


6.32li;;?« 
5.481";?! 


1 qC+O.OS 

1.00_„ 21 

1 17+0.13 
■'-.-'■'-0.18 


cluster 


18 21 58.3 


+68 18 41.2 


3.66+0.39 


1.692+0.008 


53tl 




QSO 



First line refers to the Einstein-de-Sitter model with Hq = 50 km s ^ Mpc ^ , ^Im = 1, and f^A = 0, while second line refers 
to the cosmological concordance model with Ho = 70 km s~^ Mpc~^, Q.m = 0.3, Q.a = 0.7 



